Strong-Coupling Superconductivity in NaFei-^Co^As: the Eliashberg Theory and 

Beyond 
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We study the normal-state and superconducting properties of NaFei-^Co^As system by specific 
heat measurements. Both the normal-state Sommerfeld coefficient and superconducting conden- 
sation energy are strongly suppressed in the underdoped and heavily overdoped samples. The 
low-temperature electronic specific heat can be well fitted by either an one-gap or a two-gap BCS- 
type function for all the superconducting samples. The ratio 7atT' ( ? / H^(0) can nicely associate the 
neutron spin resonance as the bosons in the standard Eliashberg model. However, the value of 
AC/T c yN near optimal doping is larger than the maximum value the model can obtain. Our results 
suggest that the high-T c superconductivity in the Fe-based superconductors may be understood 
within the framework of boson-exchange mechanism but significant modification may be needed to 
account for the finite-temperature properties. 
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Strong-coupling superconductivity in the conventional 
superconductors can be well described within the frame- 
work of Eliashberg theory [1J, where electron Cooper 
pairs are mediated by virtual phonons or some other 
bosons. Since spin fluctuations may act as the medi- 
ating bosons for electron pairing and superconductivity 
h is important to determine if the Eliashberg-based 
theory can understand the transport and magnetic prop- 
erties of unconventional superconductors [J] . For copper 
oxides, this is difficult due to the plethora of phases com- 
peting with superconductivity and the d-wave nature of 
the superconducting gap symmetry. The Fe-based su- 
perconductors may offer a better opportunity to test the 
suitability of the Eliashberg theory due to the s-wave 
nature of the superconducting electron pairing and the 
Fermi- liquid- like normal states [B| . 

In the standard Eliashberg theory, the superconduct- 
ing electron Cooper pairs are mediated by bosons with 
an average energy of w/„. For a 6- function electron-boson 
spectral density a 2 F(w) = AS(w — we), we have wi n 
= we- The ratio of T c /wi n representing the coupling 
strength is related to two important dimensionless pa- 
rameters 7ArT 2 /ff 2 (0) and AC/T C ~/ N , where j N , H c (0) 
and AC/T C are the normal-state Sommerfeld coefficient, 
the thermodynamic critical field at zero temperature, and 
the specific heat jump across T c , respectively For con- 
ventional superconductors, these two ratios can be solved 
analytically through 
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should be simultaneously satisfied for a given supercon- 
ductor. 

Recently, the bosonic spectrum is found in the tun- 
neling measurements on the Fe-based superconductors 
[^-Q. It is thus important to determine to what extent 
the standard Eliashberg theory holds by checking the va- 
lidity of Eqs. (1) and (2). The thermodynamic proper- 
ties of the Fe-based superconductors have been measured 
in many systems [9|-|20(. In Ba 6 K 4 Fe 2 As2, a calcula- 
tion based on the Eliashberg model considering multiple 
bands is able to quantitatively reproduce the experimen- 
tal results based on the assumption that spin fluctuations 
are electron pairing mediating bosons [12J. Recently, a 
very sharp neutron spin resonance is found in supercon- 
ducting NaFeo.955Coo.o45As [5l|]. The mode, centering 
at the in-plane antiferromagnetic wave vector, is strictly 
two-dimensional in the reciprocal space, which leads to 
an easy way of considering a 2 F(w) and hence wi n . There- 
fore, the NaFei^Co^As system may be suitable to check 
the Eliashberg theory. 

In this paper, we report a comprehensive study on 
the electron-doping evolution of the specific heat in 
NaFei-zCozAs. The measured value of jnT 2 /H 2 (0) is 
consistent with that obtained from Eq. (1) by assuming 
that wi n is equal to the neutron spin resonance energy 
2l| . However, the value of AC/T c ^n reaches up to 3.7 



We see that these two ratios has a linear relationship and 



near optimal doping, which is much larger than the max- 
imum value of Eq. (2). Our results suggest that the 
high-T c superconductivity in the Fe-based superconduc- 
tors may be understood within the conventional boson- 
exchange mechanism but the finite-temperature proper- 
ties should be revised around the optimal doping. 

Single crystals of NaFei-^Co^As were grown by the 
self-flux method as reported previously [22|. The sam- 
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FIG. 1: (Color online) (a) Phase diagram of NaFei-xCozAs 
that shows T c (black solid square), TV (red open circle) and 
T s (blue open triangle). The dash lines are guided to the eye. 
(b) Specific heats of some samples plotted as C/T vs T. (c) 
Fitted results of several models on the specific heat of the x = 
0.2 sample. The differences between the data and each model 
are given in (d). 



pies were attached onto the heat capacity pucks in the 
glovebox and transported within a sealed bottle to avoid 
the sample quality change 22j. The time that the sam- 
ples were exposed to air during the installation of the 
puck was less than 1 minute. The specific heat was mea- 
sured by the PPMS from Quantum Design. 

The phase diagram of NaFei-zCozAs is very similar 
to other iron pnictides with a long-range AF order in 
the parent compound and a dome-like superconducting 



regime as shown in Fig. HJa). The struc- 

tural transition temperature T s and magnetic transition 
temperature XV determined from the resistivity measure- 
ment are similar to those reported in the other literatures 
[IS 

|24{. The T c is obtained from the specific heat mea- 
surement and it is set to zero for those that exhibit no 
superconducting jump despite the fact that the resistiv- 
ity goes to zero in some samples [2tSj . Therefore, the 
superconducting dome plot in Fig. Ufa) only includes 
the samples that show bulk superconductivity. 

The raw data of specific heat are plotted in Fig. 
QJb) for some of the samples. It is clear that the 
phonon contribution varies a lot for the different Co 
doping, which makes it impossible to use the spe- 
cific heat of non-superconducting samples (e.g, x=0 or 
0.2) as a reference to determine the electronic spe- 
cific heat of superconducting samples as done in some 



other materials [13, [lj, l27|. To understand the spe- 



cific heat of the non-superconducting samples, we con- 
sider a Debye plus Einstein model assuming the to- 
tal specific heat to be C = jnT + Cd + C E , where 
C D = A D (T/T D ) 3 J Td/T x 4 e x /(e x - l) 2 dx and C E = 
A E (T E /T) 2 e TE/T /(e TE/T -l) 2 are the specific heats from 




FIG. 2: (Color online) The electronic specific heat (black open 
circle) of serial samples obtained as described in the text. The 
red lines are the fitted results of one-gap BCS function except 
for the x=0 sample which shows no superconducting jump. 
The blue line in (h) is the subtracted data between and 9 
Tesla. 



the Debye and Einstein models, respectively. Fig. [TJc) 
shows the fitting results on the x=0.2 sample for the De- 
bye+Einstein model and some other models. Fig. []Jd) 
further gives the differences between the raw data and 
the fitting results of various models, which unambigu- 
ously shows that the Debye+Einstein model gives the 
best fit to the data. We note that the parameters such 
as the Debye temperature and Einstein temperature in 
the fitting may not reflect the real phonon physics in this 
system. Since it is only possible to fit the normal-state 
and low-temperature data of the superconducting sam- 
ples ( assuming that the superconducting gaps are fully 
opened), we also test the above method by removing the 
x = 0.2 data between 3 K and 20 K in the fitting and the 
result is consistent with that fitted with the whole tem- 
perature range. The value of jn is manually adjusted for 
the superconducting samples to make sure that the en- 
tropy is conserved. In addition, the residual Sommerfeld 
coefficient 70 is obtained by fitting the low-temperature 
specific heat with C = 70 T + f3T 3 . 

The subtracted electronic specific heats of the super- 
conducting samples are shown in Fig. [2] All the data 
except for the x = 0.015 can be well fitted by the two-gap 
BCS expression of the specific heat ( C = A\Cbcs{^i) + 
A2Cbcs{&-2)) as shown by the solid lines [131] . Only one 
gap is needed to fit the x = 0.015 data. It should be 
pointed out that we cannot rule out the existence of nodes 
or highly anisotropic gaps 2^, 29[ due to the limitation 
of our model. Fig. shows the doping dependence 
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FIG. 3: (Color online) Doping dependence of (a) 7at ( black 
solid squares ) and 70 ( red open circles ), (b) the condensation 
energy, (c) 2A(0)/fc s T c and (d) Ai/(Ai +A2). All the dashed 
lines are guided to the eye. 



of 7at and 70. Contrary to that in Ba(Fei_ a; Co a ;)2As2 
27], 70 is much smaller than 77V for all the supercon- 
ducting samples, suggesting that most of the electrons in 
NaFei-^Co^As are condensed at K. With increasing 
Co doping, 7tv quickly increases and reaches its peak at 
the optimal doping with x = 0.025. Such suppression in 
the underdoped regime is most likely due to the opening 
of the SDW gap 27J. Further increasing Co above 0.05 
rapidly reduce to a very low value for heavily over- 
doped samples ( 7^ = 2.3 mJ/mol/K 2 for the x = 0.08 
sample ). Surprisingly, the goes back to more than 3 
mJ/mol/K 2 for the x > 0.1 samples. 

The doping dependence of may be strongly associ- 
ated with the SDW gap and pseudogap as observed by 
the STM 30]. Since Co doping only shifts the Fermi 
level without significantly changing the band structures 
[3~0 ] , we can quantitatively estimate the effect of the two 



gaps. Taking that N(0)(l+A)=0.42Wn [17| where n = 
3 and N(0) « 0.53 states/eV/atom [32j, we get that the 
coupling parameter A is about 0.9 for the x=0.025 sam- 
ple. Such value is reasonable considering that no pseu- 
dogap is found near the optimal doping [30]. Suppos- 
ing A does not change below 0.1 Co doping, we estimate 
that the DOS is suppressed about 60% for both the x 
=0 and x = 0.06 samples, which is consistent with the 
STM results [13]. We note that the suppression of 
in overdoped NaFei-zCo^As is much larger than that 



in Ba(Fei_ a; Co a ;)2As2 [27j, which suggest that the latter 
may have a different origin. Increasing Co doping above 
0.1 results in the depinning of the large "V" -shaped fea- 
ture and thus the disappearance of the pseudogap [30I ]. 
Our measurements on the 0.1 and 0.2 samples show that 
A is close to zero assuming that there is no suppression of 
DOS at Fermi level, which accords with the fact that the 
system is close to a normal metal with weakly coupled 
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FIG. 4: (Color online) (a) Theoretical result of (7jv — 
7o)Tc /H% (0) calculated from Eq. (2) as shown by the red 
line. The values of conventional superconductors fall into the 
shaded area. The open circles represent the values obtained 
in this paper which give the corresponding T c /ui n (b) The 
T c dependence of the resonance energy En and uji„. The 
solid line is guided to the eye. (c) The doping dependence of 
AC/T c ( 7iv - 70). (d) The corresponding ( 7JV - 7o)T 2 /tf 2 (0) 
and AC'/Tc(-yN — 70) for each sample. The solid line is the lin- 
ear relationship between these two values as calculated from 
Eq. (1) and (2). 



electrons in this doping regime. 

Fig. E^b) plots the doping dependence of the con- 
densation energy E c , which is obtained through E c = 
— Jq" C/TdTdT. It is clear that E c is much smaller 
at either the x = 0.015 or x = 0.06 sample, consistent 
with the fact that the DOS at the Fermi level is strongly 
suppres sed due to either the SDW gap or the pseudogap 

Mm- 

The doping dependences of the two superconducting 
gaps and the relevant ratio of the small gap are shown 
in Fig. [S^c) and EK<4) respectively. The values of the 
larger gap above x = 0.025 are more or less consistent 
with the results of ARPES and STM experiments where 
only one gap is observed 1341 [3BI ] . The existence of 
the smaller gap and its increasing contribution to the 
electronic specific heat are missing in those experiments. 
Since the tunneling matrix element of the M-centered 
bands may be strongly suppressed for a good surface in 
the STM experiment [361 ] , the small gap may exist around 
the M point with a non-zero k z value [jyj]. For the x = 
0.015 sample where the AF order presents, it is not clear 
why a much smaller gap is obtained from the specific heat 
data dHHi]. 

The two dimensionless ratios in Eq. (1) and (2) can 
be derived from the above experimental data. To elimi- 
nate the effect of the residual electronic specific heat, we 
replace 7 at in the ratios to 7 at — 70. Fig. 4(a) shows 
the doping dependence of (7^ — 7o)T c 2 /i? 2 (0) by taking 
H c (0) 2 — 8irE c . For the conventional superconductors, 
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the ratios of many materials are within the shaded area 
centering the red line in Fig. |4|a) calculated by Eq. (1) 
[l[ . By assuming that it is on the other side of the curve 
in our case (a value of T c /uji n much smaller than 0.2 will 
result in a bosonic energy that has not been observed 
in other experiments), we are able to obtain T c /w/„ ( 
the open circles in Fig. |4ja) ). Fig. 4(b) shows the T c 
dependence of uji n , which gives uji n = 3.38/cbT c . Inter- 
estingly, the resonance energy in the NaFeo.gssCoo.tHsAs 
is very close to the value of uji n [2l|, suggesting the res- 



onance mode may play as bosons in the superconductiv- 
ity of NaFei-zCozAs. It will be interesting to compare 
the resonance energy in heavily underdoped and over- 
doped samples with the w;„ obtained here. The large 
values of T c /uii n suggest the strong-coupling nature of 
NaFei-zCozAs noting that the ratio only extends up to 
about 0.25 in the conventional superconductors. We note 
that the above uji n can also give a value of 2A(0)/fc_B?c 
larger than 5 as suggested by the equation (4.1) in Ref. 

a. 

The doping dependence of the normalized specific heat 
jump AC/Tc^jy — 70) is shown in Fig. SJc). In the case 
where the two-gap BCS function cannot give a good fit 
near T Cl a sing-gap BCS function is used to just fit the 
data near T c to obtain an accurate specific heat jump. 
A dome-like feature is seen as that in Ba(Fei_ :r Co :c )2As2 
27]. Surprisingly, the largest value of AC/T c (jn — 7o) 
is 3.7 for x = 0.025, which is much larger than those 



found in other systems 
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While the values of 



(7jv — lo)Tc / H%(0) seem to be reasonable, Eq. (2) fails 
to calculate AC /T c {^n — To) at the optimal doping. This 
is clearer by plotting these two ratios together as shown 
in Fig. SJd). The data falling on the solid line suggest 
that they can be calculated from each other according to 
Eq. (1) and (2). It is clear that strong deviation occurs 
near the optimal doping. Since (jn — 7o)r c 2 /-ff^(0) is 
close to each other except for that of x = 0.06, such de- 
viation is not due to the insufficiency in calculating the 
two ratios by the perturbation method [ij. For the x = 
0.06 sample, it is possible that a more accurate method 
may give a better result or the pseudogap-like phase may 
has something to do with the mismatch. 

A large specific heat jump in the Eliashberg model is 
a result that the superconducting gap opens up more 
rapidly just below T c than it does in the BCS theory 
[l[, which will only give a maximum value of about 3 as 
seen in Eq. (2). In the case of NaFei-^Co^As, one may 
has to consider a very weak temperature dependence of 
the gap [HI]. On the other hand, the ratio -f N T^/H^(0) 
is associated with the condensation of the Cooper pairs 
at zero K, which may not contradict with what happens 
near T c . After all, a very important assumption in the 
strong-coupling theory is that the boson spectrum is fixed 
while the spin fluctuations in the Fe-based superconduc- 
tors strongly evolve with changing temperature. 

It is also found in the heavy-fermion materials CeCoIns 



37| and CeIrSi 3 [38| that AC/T c j N is lager than 4. This 
is consistent with theories associated with the strong lo- 
calized spin fluctuations (39l - l4lj | , indicating that spin fluc- 
tuations may result in the largest enhancement of the 
specific heat jump near optimal doping 42j. 

In conclusion, we test the validity of the Eliashberg 
formalism in the NaFei-^Co^As system by deriving 
7ArT c 2 /iJ c 2 (0) and AC/T c (j N ~ 7o ) from the specific heat 
measurements. Our results show that while the former 
value is nicely associated with the neutron spin resonance 
through Eq. (1), the latter value is beyond the Eliash- 
berg theory near optimal doping. Therefore, the pairing 
mechanism in NaFei^^Co^As may be understood within 
the boson-exchange mechanism but the disappearance of 
the superconductivity near the optimal doping should be 
considered with significant modification of the theory. 
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